Coherence in a coupled network: Implication for brain function 
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In many body systems, constituents interact with each other, forming a recursive pattern of mutual 
interaction and giving rise to many interesting phenomena. Based upon concepts of the modern 
many body theory, a model for a generic many body system is developed. A novel approach is used 
to investigate the general features in such a system. An interesting phase transition in the system 
is found. Possible link to brain dynamics is discussed. It is shown how some of the basic brain 
processes, such as learning and memory, find therein a natural explanation. 
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A major task in physics research is to search for gen- 
eral principles behind various nature processes. One of 
the greatest advances in the last century is understanding 
of how a many body system manifests and behaves, lead- 
ing to many important phenomena including the multiple 
scattering of both quantum and classical waves, the quan- 
tum Hall effect, the transition between super-conducting 
and normal metal states and so on. It has become a daily 
experience that macroscopic objects arise from interac- 
tion between individuals in a many body system which 
may be of microscopic scale. 

In this way, we were led to a fundamental question: are 
there general features which can be said about the man- 
ifestation of many body interactions. Addressing this 
question, many profound concepts have been developed 
and matured over the past half century. To name a few, 
■this includes the phase transition, invariance or symme- 
,try, and symmetry breaking; these concepts fabricate the 
■quilt of the modern many body theory |l[. An intricate 
, inter-connection among these is stated by the Goldstone 
■theorem. A symmetry breaking will result in not only a 
'phase transition, but also the appearance of some sort of 
long range order, or synchronization, or coherence, de- 
pending on what subject is under scrutiny. For example, 
the phase transition between the liquid and solid states 
jis a result of breaking the continuous space invariance. 
■Once the symmetry breaks, phonon modes appear to re- 
I fleet the collective behavior of the system. Although they 
■can be systematically obtained in the realm of the quan- 
,tum field theory, for instance, from the Ward-Takahashi 
■ relation, these phonon modes are not of quantum origin. 
This observation leads us to believe that the long-range 
macroscopic coherence or order does not necessarily re- 
quire a quantum cause. 

Being a most complicated classical hierarchy kingdom, 
human brain functions in such a way that it apparently 
involves many interacting elements. Natural questions 
may arise as whether the concepts of the many body 
theory can be applied to the problem of brain, and what 
general features can be depicted about the brain function. 
Physical research on brain function is not merely out of 
the scientific curiosity, but rather the problem is one of 
the most basic research subjects and has long been an 
outstanding inquiry from various scientific communities. 



In view of the complexity of a brain, any task on brain 
function seems unfeasible at first sight. Towards this, we 
may draw attention to the study of gas systems. Even in 
a simplest perfect gas system, it is absolutely unrealistic 
to trace the trajectory of any individual gas molecules in 
the realm of Newton's mechanics. Such a seemingly in- 
tractable problem, however, has been studied successfully 
in terms of statistical mechanics. Therefore, whether a 
most difficult problem can be formulated and has cer- 
tain answers crucially depends on from which angle and 
at what level the problem is investigated. Searching for 
a general description of nature phenomena without hav- 
ing to tangle in the mist of complexity and to know the 
details of each consisting component in complicated sys- 
tems is in fact one of the major tasks for physicists. 

Motivated by our recent work on wave localization |^ , 
here we formulate a general model for interacting many 
body systems. It will be clear that this model can rep- 
resent many important physical situations, and it may 
also be applicable to model brain systems. The moral is 
that as long as there is certain consensus about the brain 
function, a general underlying principle can be identi- 
fied. Then a new approach will be proposed to study the 
model. Numerical examples will be presented, then fol- 
lowed by a discussion on their implications for brain func- 
tion. It will become clear that some of the basic brain 
processes, such as learning and memory, find therein a 
natural explanation. Undeniably, the present investiga- 
tion has to be primitive due to limited information and 
the type of information available about the brain. Nev- 
ertheless, the study does provide a most general insight 
to the problem of brain function with respect to features 
involving many body interaction. 

Without compromising generality and for the sake of 
simplicity, we start with considering a system consisting 
of N electric dipoles, in the hope of setting up a generic 
model for many body systems. These dipoles can be ei- 
ther randomly or regularly distributed in a space which 
can be either one, two, or three dimensional. The former 
forms random arrays of objects, while the latter would 
form a crystal lattice allowing band structures to appear. 
The dipoles interact with each other through transmit- 
ting electrical waves when each dipole oscillates. 

The equation of motion for any dipole can be easily de- 



duced from the electro- magnetic theory Neglecting 
terms higher than the second order, for the i-ih dipole, 
for example, the equation of motion for the dipole oscil- 



lation is 
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In these, Qi is the oscillation displacement with uji being 
the natural frequency and 7i being the damping factor, 
Zi is the unit vector along the dipole direction of the i-th 
dipole, fij is the unit vector pointing from the j-ih. dipole 
to the z-th dipole, the couphng constant Cy — ■ 
In the present paper we ignore the rotational motion of 
the dipoles, which can be added with ease. 

The second term on the RHS of Eq. (|l|) denotes the 
nonlinear interaction. For the first step analysis, this 
effect may be dropped. Equation can formally be 
generalized as 



(2) 



where "01 is the physical quantity under consideration; 
in the dipole case ■i/'i is the displacement, Fi refers to 
the external stimuli at the i-th entity, is an operator 
transferring interaction between the i-th and j compo- 
nent which usually can be written in terms of a Bosonic 
propagator. The operator C{d) is called divisor obtained 
when the many body interaction is not present So 
far, in nature there are two known types of divisors: 



Type I: C{d) 
Type II: C{d) 
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Type I usually refers to the boson like systems, while 
Type II divisor usually refers to the fermion like systems. 
In the above w(V) and e(V) refer to the operators that 
give the energy spectra in absence of interactions, i. e. 
the bare energy. When dissipation is present, a term like 
7^ should be added to the divisors. 

Equation (|^) accommodates a variety of problems of 
great interest. To name but a few, the vibration of the 
lattice formed by atoms is described and multiple scat- 
tering of waves in a system consisting of many scatterers 
are described by a version of Type I equation. For an 
electronic system, in terms of the Local Combination of 
Atomic Orbits, the behavior of the system is governed by 
a Type II equation. Other examples can be referred to 
in fl. 



In short, Eq. (g) represents a wide range of many body 
problems. It reflects a basic many body interaction pic- 
ture. The interaction among the constituents is mediated 
by waves. This is common in nature. We would like to 
explore some new features embedded in Eq. (||). We con- 
sider wave propagation in a Type I system. The incident 
wave will be interfered by radiated waves from each en- 
tity. The energy flow is J ~ Re[M*(— iV)?/]. Writing the 
field as Ae'^, the current becomes A^VO, a version of the 
Meissner equation. This hints at that when 6 is constant 
while A 7^ 0, the flow stops and energy will be localized 
or stored in space. This phase transition is condensation 
of modes in the real space, in contrast to the Cooper 
pair condensation in the momentum space in the super- 
conductivity. This is useful as it implies that energy can 
be stored in certain spatial domians; as will be clear, it 
links to some basic features of brain function. 

We solve Eq. (||) for Type I situations in the frequency 
domain. There are N bodies in the system, and they 
are excited by a point source located in the middle of 
these bodies. To further simplify our discussion, we as- 
sume that all elements in the system are identical and 
the interaction is isotropic; the study of anisotropic cases 
is possible but comuptationally expensive and some ex- 
amples have been documented elsewhere |^]. Explicitly, 
from Eq. (^) the governing equation can be written as 
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where C is the coupling constant, and c is the wave speed. 
Strictly, for 2D {d = 2), the propagator should be in a 
form of the zero-th order Hankel function of the first kind 



We write 



(4) 



where we dropped the common time factor e~*"*. The 
total wave in the space is the summation of the direct 
wave from the source and the radiated waves from all 
the bodies. For each phase Oi, we define a phase vector 
such that 



cos tirCx 
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This phase vector can be placed on the respective entity 
and put on a 2D plane spanned by etr, e^. The amplitude 
Ai refers to the energy distributed among the entities. 

Numerical computation has been carried out for ID, 
2D and 3D situations respectively. The results are sum- 
marized as follows. (1) In the 2 and 3 dimensional cases, 
for sufficient large number density of bodies, e. g. exceed- 
ing roughly in 3D, and for small damping factor, 
the energy is trapped and stored inside the system for a 
range of frequencies located above the nature frequency. 



More explicitly, at a frequency in that range, the trans- 
mitted wave from the exciting source is localized inside 
the system. There is no energy flow, as supported by 
the appearance of a global coherence in the phase vec- 
tors defined above, i. e. all phase vectors point to the 
same direction and phase is constant. (2) When the 
damping increases, transferring more and more energy 
to other forms, the coherence phenomenon starts to dis- 
appear gradually. (3) Decreasing the density of the con- 
stituents also allows the energy to be released from the 
system, meanwhile the coherence disappear. Another ef- 
fect of decreasing density is to narrow the trapping fre- 
quency range. (4) Outside the trapping frequency range, 
no energy can be stored in the system. Therefore ad- 
justing the frequency will crucially tune whether the sys- 
tem can store energy. (5) For ID cases, energy is stored 
near the source for all frequencies with any given amount 
of random placement of the entities. (6) Last, perhaps 
most importantly, all these features are not sensitive to 
any particular distribution of the constituents. In other 
words, these features remain unchanged as the configu- 
ration of the system is varied. 

These features are illustrated by a 2D situation shown 
in Fig. 1 Here is shown that no energy is stored for 
///o = 0.98039 and 5.8824, but the energy becomes 
trapped for ///o = 1.9608. Clearly there are phase 
transitions from ///o = 0.98039 to 1.9608, and from 
1.9608 to 5.8824. We stress that these features do not 
require quantum origin. Here /o refers to the nature fre- 
quency. It is also clear from Fig. 1 that when the energy 
is stored, a striking coherence behavior appears in the 
system. Similar features are also observed for one and 
three-dimensional situations. 

In order to characterize the phase transition, we define 
an order parameter as A = X^i^i '^i ■ The magnitude 
of the order parameter represents the degree of coherence 
in the system. Obviously, when all the phase vectors 
are pointing to the same direction, the order parameter 
will equal one. The order parameter is plotted for one 
random realization as a function of the reduced frequency 
///o in Fig. 2. The two phase transitions appear around 
///o = 1.6 and 3.6. The apparent fluctuations below 
///o — 1.5 and above 4 are due to effects of random 
distribution and flnite sample size. 

Now we are ready to discuss the linkage to brain func- 
tion. Although brain may be the most complex system 
in the nature, still several general observations (dogmas) 
can be stated. First, brain is a system of many body. 
Any single neuron should not be signiflcant for the whole 
brain, but rather the patterns of activity of clusters are 
important. As early as in 1967, Ricclardi and Umezawa 
first proposed to explain the brain dynamics in the con- 
text of many body theory |^ . Anderson |0] also discussed 
some concepts in the many body theory with connection 
to living matter. Their work has directly or indirectly 
stimulated and continues to inspire much research along 
the line (e. g. [p|-pT| ) . A special area in brain research has 
emerged and is known as quantum brain dynamics ]l^ . 



Second, experiments have indicated that there are coher- 
ence phenomena in brain function. Neurons or neuron 
clusters at different spatial points act collectively, per- 
forming a specific task. Such coherence is also observed 
for microtubles in nerve cells ||l^ and other biological 
systems |Q . To effectively performing long range collec- 
tive tasks requires the communication to be mediated by 
waves and the brain system to be boson-like. Third, brain 
not only accumulates but stores energy, corresponding to 
learning and memory processes respectively. As memo- 
ries in a healthy brain are flashed in a regular pattern, 
it is reasonable to believe that the memory as a form of 
energy is stored in a coherent manner. Frolich was the 
first to link energy storage with long range coherence in 
certain biological structures jlj] . Fourth, the brain func- 
tion can resist reasonable damages which may cause the 
brain structure to dislocate to a rational extent. This 
implies that the brain function can be executed in rel- 
atively disordered environments. Lastly, brain can age, 
losing memory and becoming dying. 

All these important features of brain can be well ex- 
plained in the realm of the above generic model. Al- 
though it is not possible to discern the exact entities 
that carry out brain function at this stage, a few pos- 
sibilities have been mentioned in the literature, such as 
neurons or corticons, proteins in microtubules and so on. 
Even water molecules are also mentioned as a core role 
in memorizing and storing [p^ . Putting details aside, 
what is important is only a quantity somehow related to 
the activity of brain. Actually, it is known that neu- 
rons do not seem to be the only fundamental units of 
the brain; or, at least, together with them some other 
elements, as for instance glia celles, may play an impor- 
tant role |l^. No matter what they exactly are, the 
necessity for stability requires the constituents to follow 
Type I equation. The communication between entities 
should be mediated by waves to ensure effective infor- 
mation exchange. Based on the above results, we argue 
that the brain consists of many subsets of Type I enti- 
ties, each subset has a different natural frequency, and 
accomplishes different duties. Different subsets will store 
specific information depending on the spectra of external 
stimuli. Once the stimuli fits into a trapping range of 
a specific subset, the energy of characters of the stim- 
uli will be stored or imprinted into that subset in a co- 
herent manner, as indicated by Fig. 1. In other words, 
different subsets response to respective stimuli and store 
energy accordingly, corresponding to learning and mem- 
ory processes. When the density of entities is not high 
enough or there is no enough active entities, however, 
the energy cannot be stored for long time. This reveals 
as the situation of short memory; this may be seen, for 
example, in the early development of children's brains. 
Long term memory requires a sufficient amount of ac- 
tive entities and reasonable supply of stimuli. Moreover, 
as mentioned above, the feature of long range coherence 
and energy storage remains unchanged as the entity con- 
figuration varies; here we exclude the extreme case that 



all entities collapse into each other. This would explain 
the brain's resistance to reasonable damages. Damages 
may reduce the active entities and alter the structure of 
the system. Either varying the number of active entities 
or the natural frequency of entities or both produces the 
different states of the brain and controls the energy re- 
lease. Afcmorics may also be lost and disoriented when 
the the damping effect becomes prominent. The possible 
contributions to the damping effect include those from 
aging effect and diseases. 

In summary, while there are many unanswered ques- 
tions pertinent to brain function, necessitated by lim- 
ited knowledge about brain, it is almost certain that the 
brain is a many body system which can accumulate and 
store energy. When needed, the stored energy can be 
released to accomplish missions. The communication be- 
tween different parts in brain is likely mediated by certain 
brain waves, as no other means can convey information 
more effectively than waves. As long as these general 
discernments hold, given its generality the present model 
expects to play a role in guiding or stimulating further 
experimental and theoretical studies. 
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FIGURE CAPTION 

Figure 1 Energy distribution and diagrams for the 
phase vectors in a 2D random configuration. Right: 
Energy distribution; the geometrical factor has 
been scaled out. Left: the phase diagrams for the 
phase vectors. 

Figure 2 Order parameter versus ///o- 
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